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1 . Supplementary Methods 

1.1. Hi-C libraries, reads mapping and filtering and generation of contact 
matrices. 

Hi-C libraries from T47D cells treated or not with the Pg analogue R5020 for 60 min 
were generated according the previously published Hi-C protocol with minor 
adaptations (Lieberman-Aiden et al. 2009). Hi-C libraries were performed in both 
conditions using Hindlll and Ncol restriction enzymes to generate two independent 
biological and technical replicates. Briefly, cells were fixed with 1% formaldehyde 
during 10 min at room temperature. Cross-linking reactions were stopped by addition 
of glycine (0.125 M final). Cells were scraped and nuclei were prepared as described 
previously. Chromatin digestion, labelling and ligation steps were performed 
according the original protocol (Lieberman-Aiden et al. 2009). After deproteinisation, 
removal of biotinylated free-ends and DNA purification, Hi-C libraries were controlled 
for quality and sequenced on an lllumina Hiseq2000 sequencer. Paired-end reads 
were processed by first aligning to the reference human genome (GRCh37/hg19) 
using BWA. Reads filtered from analysis include those that were not uniquely 
mapped, mapped more than 500 bp from relevant restriction sites, had bad sequence 
quality {e.g. 40% or more of the bases with Sanger PHRED quality < 2) or bad BWA 
mapping quality (< 20), or were located in regions classified to have exceptionally 
high sequence depth (top 0.1%) by the 1000 Genomes project's data. To offset bias 
introduced by PCR amplification of the sequencing library, only one of the duplicated 
pairs were used for subsequent analyses. Datasets normalized for experimental 
biases and sequencing depth (Lieberman-Aiden et al. 2009) were used to generate 
contact matrices at 20, 40 and 100 Kb as well as at 1Mb resolutions. The 
Supplementary table 1 summarizes the number of interactions obtained for each 
dataset. 

1.2. Genome segmentation in TADs. 

The genome was segmented in TADs by the TADbit program that includes a 
change-point algorithm for the detection of TAD borders inspired by methods used to 
detect copy number variations in CGH experiments (Pique-Regi et al. 2008). Briefly, 
the optimal segmentation of the chromosome in k TADs is computed by maximum 
likelihood for every k, after which a Bayesian Information Criterion selects the best 
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model. The hypotheses underlying the model are that the number of interactions 
between two loci has a Poisson distribution, of which the average decreases as a 
power law function of their separation (in base pairs). Each TAD corresponds to a 
vertical slice of the Hi-C interaction matrix, where the parameters of the model 
mentioned above are constant. The TADbit program uses an internal normalization 
and thus runs on raw interactions. The parameters are corrected for local biases 
(most notably G+C content, availability of the restriction enzyme site and repeat 
coverage) (Lieberman-Aiden et al. 2009) that can affect the read count of the bin on 
column / and row j by taking the ratio of the sum of reads in column / multiplied by 
the sum of reads in row j. The detail of the implementation of TADbit will be further 
detailed elsewhere (Serra et al., manuscript in preparation). 

1.3. Epigenetic data collection and analyses. 

MNase-seq, DNase l-seq and ChlP-seq experiments (H3K4me1, H3K4me3, H2A, 
H4, RNA polymerase II, progesterone receptor, H3K9me3, HP1 y) in T47D cells were 
described previously (Ballare et al. 2013; Vicent et al. 2013). Additional ChlP-seq 
experiments for CTCF, H3K36me2, H3K27me3, H3K14ac and H1.2 were performed 
in similar conditions using the following antibodies: 07-729 (Millipore), 07-369 
(Millipore) 39155 (Active Motif), 07-353 (Millipore) and ab4086 (Abeam), respectively. 
All reads were processed by aligning to the reference human genome 
(GRCh37/hg19). MNase-seq, DNase l-seq and ChlP-seq signals normalised for 
sequencing depth were summed in windows of 100 Kb. To generate plots in 
Supplementary Figures 2 and 3, bins correspond to a fifth of the TAD size (hereafter 
mentioned as sub-segments). Summed reads for these different window sizes were 
divided by the corresponding signal obtained for an input DNA of T47D cells to 
determine the normalised signal over input enrichment/depletion. Progesterone- 
induced enrichment or depletion of mark content were determined for 100 Kb bins or 
for the whole TAD as the ratio of sequencing-depth normalised read counts before 
and after treatment. 

1.4. Role of TAD borders in demarking chromatin blocks. 

To analyse the role of the TAD borders in limiting epigenetic blocks of individual 
marks, the opposite of the absolute difference of the signals, normalised as above, 
was calculated for two consecutive sub-segments (with each TAD divided into 5 
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segments of equal size) on 3 consecutive TADs. Therefore, the higher the score is, 
the more similar the signals between consecutive segments. These scores were 
computed genome-wide, and the quartiles of the 14 consecutive values were 
calculated to generate the plots in Supplementary Fig. 3a and c. To control that 
the observed differences were not due to biases in the sizes of the sub-segment, 
and to confirm the role of the TAD borders, we shuffled the TADs chromosome-wise. 

To analyse the homogeneity of epigenetic marks within TADs, we generated 
pairwise correlation matrices between the profiles of the marks described above for 
all 100 Kb windows of the genome. From this, we computed the average correlation 
between 100 Kb windows located in the same TAD. If TADs are homogeneous in 
chromatin marks combinations, this average correlation will be higher than expected 
by chance. To estimate the null distribution of this score, we shuffled the TADs 
chromosome-wise 5,000 times and applied the same procedure. The average 
correlation obtained with the initial positions of the TADs was higher than the 
average correlation when the TADs were randomized. Finally, we repeated the same 
analysis using the differential of the normalised signal after treatment with R5020 
and the signal before treatment. 

1.5. Expression levels (RNA-seq) and hormone-induced changes. 

RNA-seq experiments were performed in T47D treated or not with 10~ 8 M R5020 (Pg) 
for 1 or 6 h or with 10~ 8 M estradiol (E2) for 6 h. Paired-end reads were mapped with 
the GEM mRNA Mapping Pipeline (v1.7) (Marco-Sola et al. 2012) using the latest 
gencode annotation version (v.18) (Harrow et al. 2012). BAM alignment files were 
obtained and used to generate strand-specific genome-wide normalized profiles with 
RSeQC (Wang et al. 2012) software. Exon quantifications summarized per gene for 
expression level determination were obtained either as normalised read counts or 
reads per kilobase per milion mapped reads (RPKM) using Flux Capacitor 
(Montgomery et al. 2010). Fold changes (FC) were computed as the log 2 ratio of 
normalised reads per gene obtained after and before treatment with hormones. To 
analyse the changes of non-protein coding regions of the genome, the number of 
normalised reads were computed for chromosomal domains that do not overlap with 
any annotated protein-coding gene. 
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In order to compare the observed and expected number of genes positively or 
negatively affected by hormone treatment and to exclude potential biased responses 
depending on low/high basal expression of genes, randomizations of protein-coding 
gene positions used in Figure 2F were obtained as follow: genes were classified in 5 
classes of equivalent sizes according their expression level. Gene positions were 
then shuffled for the 5 classes allowing conserving an equivalent average expression 
of genes per TAD and therefore avoiding the effect of potential biased distributions 
frequencies of genes positively or negatively modified depending on the basal levels 
of expression. Shuffled lists were then used to calculate the percentage of genes with 
positive or negative fold change per TAD. 

1.6. Classification of TAD according transcriptional response to Pg. 

To classify TADs according their hormone response, we calculated the average ratio 
of the number of normalised RNA-seq reads obtained after and before hormone 
treatment in the RNA-seq replicates. TADs containing more than 3 protein coding 
genes were saved for further analysis and ranked according to the average ratio 
described above. The top and bottom 10% were classified as "activated" and 
"repressed" TADs, respectively. 

1.7. Fluorescent in situ hybridization. 

T47D cells grown on slides were fixed with 4% paraformaldehyde in PBS 15 min at 
room temperature. After washes with PBS and permeabilisation with 0.2% Triton X- 
100 in PBS, slides were incubated 60 min with RNase A in SSC2X. Fixed cells were 
incubated 3 min in 0.1 N HCI and 1 min in 0.01 N HCI with 0.01% pepsin. Slides were 
denatured at 70 Q C for 8 min in SSC 2x-70% formamide and incubated overnight at 
37 Q C in a humid chamber with the probes (separately denatured 10 min at 80 Q C in 
hybridization buffer). 100 ng of each probe generated by nick translation using dUTP- 
biotin, dUTP-DIG and dUTP-fluorescein (Roche Applied Biosciences) was used per 
hybridization. Detection was performed with ant-biotin-Cy5 (Rockland), anti-DIG- 
rhodamin (Roche) and anti-fluorescein-Alexa 488 (Invitrogen). Machine optimized 
stacks were acquired on a Zeiss TCS SP5 confocal microscope. After deconvolution, 
3D rendering of the hybridised probes signals were generated for each channel, and 
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the pair-wise 3D distances between the centre of mass of these signals were 
computed using the Imaris Software. 



The BACs used in this study were obtained from a 32k library: 



RP11-758G19 


chrl 


: 26742822 


- 26905278 


RP11-443P17 


chrl 


: 26877973 


- 27069024 


RP11-973A19 


chrl 


: 27772963 


- 27993702 


RP11-667P18 


chrl 


: 28384701 


- 28571122 


RP11-318E23 


chrl 


: 28980912 


- 29129694 



1.8. Integrative 3D modelling of TADs. 

Hi-C data matrices 

Hi-C experimental data resulted in interaction counts between the loci of the genomic region 
of interest (i.e. the quantitative determination of the number of times each specific 
experimental ligation product is sequenced). We then applied an internal normalization by Z- 
scoring the sequenced raw interaction count data. First, we applied a log2 transformation of 
the raw counts and then their Z-score was calculated as: 

Zscoreij=— 

where |u and a are the average and standard deviation of the interaction counts for the entire 
Hi-C matrix. Such normalization allowed us to quantify the variability within the Hi-C matrix 
as well as to identify pairs of fragments that significantly interacted above or below the 
average interaction frequency. 

TAD representation 

Hi-C datasets obtained with Hindlll and Ncol enzymes were pooled to generate 
normalized Hi-C matrices at 20 Kb resolution. Consequently, each modelled genomic 
region was represented as a set of particles, one per each 20 Kb bin. Each particle 
had a radius of 100 nm, based on the relationship of 0.01 nm per base pair (bp)(Harp 
et al. 2000). Neighbour particles were constrained to lie at an equilibrium distance 
proportional to the sum of their excluded volume. Non-neighbor particle pairs (i.e. 
particles representing non-consecutive bins along the genomic sequence) were 
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instead assigned a distance derived from the z-score matrices. The function mapping 
the z-scores onto distances was defined as the calibration curve built by interpolating 
the highest and the lowest z-score values, with a minimum proximity distance of 200 
nm (the excluded volume of two base particles) and a maximal proximity distance for 
two non-interacting fragments, respectively, as previously described (Bau and Marti- 
Renom 2012). The forces applied to the defined restraints were also set 
proportionally to the absolute value of the Hi-C z-score observed between a pair of 
fragments. 

TAD modelling 

Each particle pair was restrained by a series of harmonic oscillator centred on a 
distance derived from the experimental data. Consecutive particles (i.e. particle pairs 
—or bins— /, i+1) were considered as neighbour particles and therefore retrained at 
an equilibrium distance proportional to the sum of their excluded volume. Non- 
neighbour particles (i.e. particle pairs —or bins— /, i+2..n) were instead restrained at 
distances calculated from a function that mapped their corresponding z-scores onto 
distances. This function corresponded to a calibration curve that was built between 
the points defined by the maximum and minimum z-score values, and two empirically 
determined minimum and maximum distances, respectively. The maximum distance 
for non-interacting loci was independently optimized for each modelled region 
(Supplementary Table 2). Consecutive particles were restrained by an upper-bound 
harmonic oscillator, which ensured that two particles could not get separated beyond 
a given equilibrium distance proportional to the sum of their excluded volume. Non- 
consecutive particles were restrained by two different oscillators: (i) an harmonic 
oscillator, which ensured a pair of particles to lie at about a given equilibrium distance 
and (ii) lower-bound harmonic oscillator, which ensured that two particles could not 
get closer than a given equilibrium distance. In both cases the equilibrium distance 
was derived by a calibration curve defined by the points corresponding to the 
maximum and minimum z-score values with a minimum and a maximum distance 
(see below). The different type of oscillator applied depended on an upper (uZ) and a 
lower (17) z-score cut-off as well as the approximation distance between two non- 
interacting particles (aD), as previously described (Bau and Marti-Renom 2012). The 
values of l/Z, \Z and aD were independently optimized for each modelled region 
(Supplementary Table 2). Neighbour particles were constrained by an upper-bound 
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harmonic oscillator centered at an equilibrium distance proportional to the sum of 
their excluded volume. Non-neighbour particles with a z-score higher than the upper- 
bound cut-off (uZ) were restrained by a harmonic oscillator, while pairs of particles 
with a z-score lower than the lower-bound cut-off (ZZ) were restrained by a lower- 
bound harmonic oscillator. Since these two harmonic oscillator aim at keeping a pair 
of particles at an equilibrium distance or further apart from a minimal distance, 
respectively, pairs of non-neighbour particles that were observed to interact with z- 
scores above the uZ parameter were kept close in space, while pairs of non-neighbor 
particles that were observed to interact with z-scores below the IZ parameter were 
kept apart. The k force applied to these restraints was set to the square root of the 
absolute value of their interacting z-scores. Finally, pairs of non-neighbour particles 
for which Hi-C data were not available were restrained based on the average z-score 
of the adjacent particles. 

Model building with IMP 

Following the steps previously described (Bau and Marti-Renom 2012), we modelled 
the 3D structure of 61 genomic regions. The selected regions included a total of 209 
TADs and covered about 267 Mb of the genome (Supplementary Table 3). Once the 
system was fully represented and the restraints between the particles were set, IMP 
generated the ensemble of solutions that best represented the input data by 
simultaneously minimizing the violations of all the imposed restraints. Due to the 
large conformational space to be explored and to the population-based nature of 3C- 
based methods, the optimization of the imposed restraints resulted in different 
configurations with similar final IMP objective function. Thus, to comprehensively 
explore the conformational space of the modelled TADs, we generated a large 
number of models for each conditions (2,000 models for each) by globally minimizing 
the imposed restraints via a combination of 10,000 Monte Carlo rounds with 5 local 
steps in a molecular dynamics simulation with a standard simulated annealing 
method. At each optimization step, the conformation was randomly changed and the 
change was accepted or rejected according to the Metropolis criteria. The driving 
scoring function that was minimized during the optimization protocol consisted of the 
sum of all the individual restraint scores between all the particles representing the 
system. The entire calculation of 122 independent simulations resulting in 244,000 
different conformational solutions (2,000 for each region and hormone treatment), 
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took about 6 days on a 200 CPU cluster computer. 
Model analysis 

To further characterize the structural perturbations on TADs by the treatment with 
Pg, we performed a series of computational analysis in a set of 1,000 selected 
models with the lowest IMP objective function, which correspond to the 3D models 
that best satisfy the initial imposed restraints (Bau and Marti-Renom 2012). 
The entire set of analysis included: 

a) Structural clustering of models. 

To structurally compare two 3D models, we used pair-wise rigid-body 
superposition minimizing the root mean square deviation between the 
superposed conformations. The comparison resulted in a 1,000x1,000 
symmetric matrix of all-against-all structural comparisons storing the 
number of particles for every pair of aligned models that align within 75 nm 
distance cut-off. The comparison matrix was then used as input file to the 
Markov Cluster Algorithm (MCL) program (Enright et al. 2002), which 
generated unsupervised sets of clusters of related structures. 

b) TAD radius of gyration. 

The radius of gyration of an object (in our case a TAD) is root mean 
square distance of the objects' parts from its centre of gravity. We 
calculated for each TAD a centre of gravity and then measured the 
distance of all points in the TAD to this centre. The average value of the 
distances corresponded to the TAD radius of gyration. The larger the 
radius of gyration is, the more open the TAD is. 

c) TAD and particle accessibility. 

The particle accessibility of a model is the accessible fraction of its particle 
to a hypothetical spherical object of a given radius (in our calculations 75 
nm). To obtain such fraction, we first build a mesh of points surrounding 
each particle in the model, which indicates its theoretical occupancy. Then, 
the spherical object is placed on all possible positions in contact with the 
model and the fraction of the mesh around the particle that is can be 
occupied by the object is considered accessible. The accessibility change 
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for a given TAD was calculated as the ratio of the sums of all particle 

accessibilities in the TAD before and after Pg. 
d) Ensemble visualization 

The UCSF Chimera package (Yang et al. 2012), a highly extensible 
program for interactive visualization of molecular structures, was used to 
produce all images of the modelled TADs. 
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2. Supplementary Tables 



Supplementary Table 1 : Summary statistics Hi-C datasets 





Hindlll -Pg 


Hindlll +Pg 


Ncol -Pg 


Ncol +Pg 


Sequenced Fragments 
Pairs 


189563259 


146045810 


186633896 


232632568 


Interacting Fragments 
Pairs* 


74015388 


48530773 


100262904 


120640304 


"Intra-chromosomal" 


29552177 


19519911 


49163452 


57721468 


"Inter-chromosomal" 


4446321 1 


29010862 


51099452 


62918836 



*: after mapping and filtering 
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Supplementary Table 2: IMP optimal modelling parameters 



Region 


Chrom 


Start 


End 


uZ 


IZ 


aD 


CC 


bin 


Bin 


1 


1 


10500 


10660 


0.0 


-0.6 


400 


0.78 


2 


1 


11315 


11590 


0.0 


-0.3 


400 


0.75 


3 


1 


1815 


1925 


0.6 


0.0 


500 


0.77 


4 


1 


3015 


3200 


0.0 


-0.3 


400 


0.68 


5 


1 


5830 


6000 


0.3 


-0.6 


400 


0.75 


6 


1 


8375 


8575 


0.0 


-0.3 


400 


0.74 


7 


1 


9070 


9315 


0.6 


-0.3 


500 


0.76 


8 


2 


11025 


11240 


0.0 


-0.3 


400 


0.68 


9 


2 


2105 


2290 


0.3 


-0.3 


400 


0.76 


10 


2 


3480 


4010 


0.9 


-0.3 


500 


0.70 


11 


2 


480 


660 


0.0 


0.0 


400 


0.78 


12 


2 


7910 


8045 


0.6 


-0.3 


500 


0.71 


13 


2 


9940 


10130 


0.0 


-0.3 


500 


0.70 


14 


3 


7515 


7920 


0.9 


-0.6 


500 


0.71 


15 


3 


8395 


8500 


0.6 


0.0 


500 


0.80 


16 


4 


1970 


2100 


0.3 


0.0 


500 


0.75 


17 


4 


355 


490 


0.0 


-0.3 


500 


0.82 


18 


4 


3990 


4165 


0.0 


-0.3 


400 


0.68 


19 


4 


4410 


4760 


0.0 


-0.6 


400 


0.63 


20 


4 


5245 


5480 


0.6 


-0.3 


500 


0.70 


21 


5 


2605 


2915 


0.0 


-0.3 


400 


0.73 


22 


5 


535 


810 


0.9 


-0.6 


500 


0.69 


23 


5 


7050 


7405 


0.9 


-0.3 


500 


0.73 


24 


6 


2680 


3200 


0.9 


-0.3 


500 


0.72 


25 


6 


3590 


3825 


0.6 


-0.3 


500 


0.69 


26 


6 


7560 


7665 


0.9 


0.0 


600 


0.72 


27 


6 


7885 


8055 


0.0 


-0.3 


400 


0.75 


28 


7 


2410 


2740 


0.6 


-0.6 


400 


0.69 


29 


7 


4390 


4610 


0.0 


-0.3 


500 


0.65 


30 


7 


7830 


7956 


0.3 


0.0 


500 


0.78 


31 


8 


2675 


2850 


0.0 


0.0 


500 


0.76 


32 


8 


455 


595 


0.0 


-0.3 


400 


0.75 


33 


8 


5195 


5345 


0.0 


-0.6 


400 


0.74 


34 


9 


0 


105 


0.6 


0.0 


500 


0.74 


35 


9 


5855 


6240 


0.3 


-0.6 


400 


0.70 


36 


9 


6960 


7060 


0.0 


-0.6 


500 


0.88 


37 


10 


1065 


1545 


0.6 


-0.3 


500 


0.67 | 
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38 


10 


2595 


2725 


0.9 


0.0 


500 


0.72 


39 


10 


295 


470 


0.6 


-0.3 


500 


0.76 


40 


10 


3850 


4065 


0.0 


-0.3 


400 


0.81 


41 


10 


4720 


4845 


0.3 


-0.3 


400 


0.82 


42 


11 


3715 


4120 


0.0 


-0.6 


400 


0.71 


43 


12 


1555 


1720 


0.0 


0.0 


400 


0.74 


44 


12 


4470 


4655 


0.6 


-0.6 


500 


0.73 


45 


12 


5115 


5290 


0.0 


-0.6 


400 


0.78 


46 


12 


805 


1075 


0.0 


0.0 


400 


0.66 


47 


13 


1195 


1305 


0.0 


-0.3 


500 


0.76 


48 


13 


1980 


2305 


0.0 


-0.3 


400 


0.68 


49 


13 


5175 


5560 


0.0 


0.0 


400 


0.66 


50 


14 


3090 


3240 


0.3 


0.0 


500 


0.75 


51 


14 


4445 


4595 


0.0 


0.0 


400 


0.76 


52 


15 


4990 


5126 


0.0 


-0.6 


400 


0.78 


53 


16 


3865 


4120 


0.0 


-0.6 


400 


0.73 


54 


18 


770 


1190 


0.9 


-0.6 


500 


0.72 


55 


20 


1570 


1670 


0.0 


-0.3 


400 


0.81 


56 


21 


1405 


1625 


0.0 


-0.6 


400 


0.71 


57 


21 


1745 


1920 


0.6 


-0.3 


500 


0.77 


58 


21 


2030 


2240 


0.6 


-0.6 


400 


0.78 


59 


23 


6470 


6575 


0.0 


0.0 


500 


0.71 


60 


23 


6860 


6990 


0.9 


0.0 


600 


0.62 


61 


23 


790 


935 


0.6 


-0.3 


500 


0.76 



CC corresponds to the correlation coefficient of a contact map based 
on the 3D models and the input HiC data for the region. High 
correlation coefficients are indicative of the bona fide representation of 
the HiC data by the models. 
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3: Summary of modelled TADs coordinates 



mm 

1 


27 
28 


1 


36,300,000 
36,900,000 
JD.UUU.UUU 


36,900,000 
38.000,000 
Jo.*>LaJ,UUu 


1,815 
1,845 


1,900 


1 

31 


31 


55 


600, D0O 
1 ,1 00,000 
Ann nrvi 








86 


2 


40 
41 
42 


1 


60,300.000 
62,400,000 
63,300,000 
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4. Supplementary Figures and legends 

Supplementary Figure 1 : T47D genome is organised in TADs. 

(A) The linear correlation between number of Hi-C contacts and number of restriction 
sites is lost due to changes in copy number allowing determining the number of 
copies of each megabase domains present in T47D genome (top left panel). A 
snapshot of contact matrix shows an example of translocation breackpoint between 
the chromosomes 7 and 15 in T47D cells (top middle panel). T47D karyotype 
deduced from the variations in copy number and translocations events determined as 
above is represented (right panel). Bottom left panel show examples of 2D-FISH on 
T47D cells metaphase spreads validating the karyotype and highlighting changes in 
ploidy, translocation events as well as normal chromosomes. (B) Distribution of 
boundaries confidence scores for TADs defined at 40 (light grey - n=3027) or 100 Kb 
(dark grey - n=2031) resolutions. (C) Distribution of TAD sizes determined at 100 Kb 
resolution (N genome=2031). 

Supplementary Figure 2: Functional characterization of TADs. 

TADs were divided into 5 sub-segments of equal size to study the distribution of 
epigenetic features according a border (b) to center (c) distribution (see 
Supplementary Methods). (A) The median percentage of Transcription Start Sites of 
protein coding genes (TSS), RNA-Polymerase II (Pol. II), and CTCF binding sites, as 
well as DNAse I hypersensitive sites (DHS) is shown. (B) The histogram represent 
the median ratio of normalised ChlP-Seq signals over input for H3K36me3, 
H3K27me3, H3K9me3 and total H4 determined in the TADs 5 sub-segments. (C) 
Normalised ChlP-Seq signals over input ratio were calculated per sub-segment for 
each TAD for the indicated chromatin components and marks mentioned. Histograms 
show, for each mark, the median normalised signal/Input observed in the 5 sub- 
segments for all TADs in the genome. 

Supplementary Figure 3: TADs are epigenetic domains which chromatin is 
coordinately modified upon Pg. 

(A) Plots show the homogeneity score of the normalised Chip-Seq signal/Input ratio 
between successive sub-segments (see Supplementary Methods) over 3 
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consecutive TADs for the chromatin marks and components listed. Lines depict the 
25 th , 50 th and 75 th percentiles (from top to bottom respectively) of the scores 
computed genome-wide. (B) Differences of +Pg/-Pg H3K9me3 Chip-Seq signal 
between consecutive sub-segments (see Fig. S3 and Supplementary information) 
over 3 consecutive TADs in the case of observed (left panel) or randomized (right 
panel) TADs borders show that transition between Pg induced changes in chromatin 
state occur preferentially at the TAD boundaries. Similar analysis was performed for 
other chromatin marks as example H3K27me3 and H3K36me2. 

Supplementary Figure 4: Homogeneous and specific response of TADs to 
steroid hormones. 

(A) The ratio of normalised reads after and before 6 hour of treatment with Pg were 
calculated for each TAD. The scatter plot shows the correlation of response to Pg per 
TAD obtained in two biological replicates (BR1 and BR2) of RNA-Seq. (B) Boxplots 
show the basal expression levels (RPKM - Log 2 ) of genes located within the three 
types of TADs. (C) Scatter plot showing the correlation of response per TAD after 6 
hours of treatment with Pg (R6) with the changes obtained after 1 hour of Pg 
induction (R1) or 6 hours of treatment with E2 (E6). 

Supplementary Figure 5: TADs respond as unit to the hormone stimulus 

(A and B) Genome browser view of RNA-Seq signal within TADs presented in Figure 
2E ((A) U469; (B) U821) highlighting the expression of non-annotated non-coding 
transcripts which correlate with the hormone induced changes in expression 
observed for the protein coding genes. 

Supplementary Figure 6: Structural changes of TADs. 

(A and B) Distributions of the changes of intra-TAD contacts proportions upon Pg 
treatment (See Extended Experimental Procedures) in the different TAD categories 
for the two independent datasets obtained with Hindi 1 1 (A) or Ncol (B). Boxplot 
whiskers correspond to 5 st and 95 th percentiles. (***)>(**)>(*) indicate P < 0.0001, 
0.001 and 0.01, respectively (Mann-Whitney test). (C) Pair-wise inter-probes 3D 
distances obtained in the models (left panel) or in situ (right panel) were plotted 
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according the genomic distances that separate them. (D) Distributions of the 
accessibility scores calculated from the models for particles containing or not a TSS. 
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Supplementary Figure 1 : T47D genome is organised in TADs. 

(A) The linear correlation between number of Hi-C contacts and number of restriction 
sites is lost due to changes in copy number allowing determining the number of 
copies of each megabase domains present in T47D genome (top left panel). A 
snapshot of contact matrix shows an example of translocation breackpoint between 
the chromosomes 7 and 15 in T47D cells (top middle panel). T47D karyotype 
deduced from the variations in copy number and translocations events determined as 
above is represented (right panel). Bottom left panel show examples of 2D-FISH on 
T47D cells metaphase spreads validating the karyotype and highlighting changes in 
ploidy, translocation events as well as normal chromosomes. (B) Distribution of 
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boundaries confidence scores for TADs defined at 40 (light grey - n=3027) or 100 Kb 
(dark grey - n=2031) resolutions. (C) Distribution of TAD sizes determined at 100 Kb 
resolution (N genome=2031). 
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Supplementary Figure 2: Functional characterization of TADs. 

TADs were divided into 5 sub-segments of equal size to study the distribution of 
epigenetic features according a border (b) to center (c) distribution (see 
Supplementary Methods). (A) The median percentage of Transcription Start Sites of 
protein coding genes (TSS), RNA-Polymerase II (Pol.ll), and CTCF binding sites, as 
well as DNAse I hypersensitive sites (DHS) is shown. (B) The histogram represent 
the median ratio of normalised ChlP-Seq signals over input for H3K36me3, 
H3K27me3, H3K9me3 and total H4 determined in the TADs 5 sub-segments. (C) 
Normalised ChlP-Seq signals over input ratio were calculated per sub-segment for 
each TAD for the indicated chromatin components and marks mentioned. Histograms 
show, for each mark, the median normalised signal/Input observed in the 5 sub- 
segments for all TADs in the genome. 
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Supplementary Figure 3: TADs are epigenetic domains which chromatin is 
coordinately modified upon Pg. 

(A) Plots show the homogeneity score of the normalised Chip-Seq signal/Input ratio 
between successive sub-segments (see Supplementary Methods) over 3 
consecutive TADs for the chromatin marks and components listed. Lines depict the 
25 th , 50 th and 75 th percentiles (from top to bottom respectively) of the scores 
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computed genome-wide. (B) Differences of +Pg/-Pg H3K9me3 Chip-Seq signal 
between consecutive sub-segments (see Fig. S3 and Supplementary information) 
over 3 consecutive TADs in the case of observed (left panel) or randomized (right 
panel) TADs borders show that transition between Pg induced changes in chromatin 
state occur preferentially at the TAD boundaries. Similar analysis was performed for 
other chromatin marks as example H3K27me3 and H3K36me2. 
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Supplementary Figure 4: Homogeneous and specific response of TADs to 
steroid hormones. 

(A) The ratio of normalised reads after and before 6 hour of treatment with Pg were 
calculated for each TAD. The scatter plot shows the correlation of response to Pg per 
TAD obtained in two biological replicates (BR1 and BR2) of RNA-Seq. (B) Boxplots 
show the basal expression levels (RPKM - Log 2 ) of genes located within the three 
types of TADs. (C) Scatter plot showing the correlation of response per TAD after 6 
hours of treatment with Pg (R6) with the changes obtained after 1 hour of Pg 
induction (R1) or 6 hours of treatment with E2 (E6). 
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Supplementary Figure 5: TADs respond as unit to the hormone stimulus 

(A and B) Genome browser view of RNA-Seq signal within TADs presented in Figure 
2E ((A) U469; (B) U821) highlighting the expression of non-annotated non-coding 
transcripts which correlate with the hormone induced changes in expression 
observed for the protein coding genes. 
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Supplementary Figure 6: Structural changes of TADs. 

(A and B) Distributions of the changes of intra-TAD contacts proportions upon Pg 
treatment (See Extended Experimental Procedures) in the different TAD categories 
for the two independent datasets obtained with Hindi 1 1 (A) or Ncol (B). Boxplot 
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whiskers correspond to 5 st and 95 tn percentiles. (***),(**),(*) indicate P < 0.0001, 
0.001 and 0.01, respectively (Mann-Whitney test). (C) Pair-wise inter-probes 3D 
distances obtained in the models (left panel) or in situ (right panel) were plotted 
according the genomic distances that separate them. (D) Distributions of the 
accessibility scores calculated from the models for particles containing or not a TSS. 
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